Metal-insulator transition on SrTiOs surface induced by ionic-bombardment 
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SrTiC>3 (STO) is one of the most popular insulating single-crystal substrates for various complex- 
oxide thin film growths, because of its good lattice match with many complex oxide films. Here, we 
show that a common thin film processing technique, argon ion-milling, creates highly conducting 
layer on the surface of STO, not only at room temperatures but also at cryogenic temperatures at 
which thermal diffusion is completely suppressed. Systematic in situ four-point conductance mea- 
surements were taken on single-crystal STO substrates inside vacuum environment. The evolution 
of metallicity out of insulating STO follows simple models based on oxygen vacancy doping effect. 
At cryogenic temperatures, ion milling created a thin - but much thicker than the argon-penetration 
depth - steady-state oxygen-vacant layer, leading to a highly-concentric metallic state. Near room 
temperatures, however, significant thermal diffusion occurred and the metallic state continuously 
diffused into the bulk, leaving only low concentraion of electron carriers on the surface. Analy- 
sis of the discrepancy between the experiments and the models also provided evidence for vacany 
clustering, which seems to occur during any vacancy formation process and affects the observed 
conductance. These observations suggest that the transport properties of films processed on STO 
substrates using energetic methods such as ion milling need to be taken with caution. On the 
other hand, if properly controlled, ionic bombardment could be used as a way to create selective 
conducting layers on the surface of STO for device applications. 

PACS numbers: 71.30.+h, 73.50.-h, 73.25.-K, 72.80.Ga 



G 

O 

i O i 

CO 

> 

CN 
CN 

in 
in 

o 



X 



I. INTRODUCTION 

Strontium titanate is one of the most studied 
perovskite-like oxides, which exhibit a variety of fea- 
tures such as high T c superconductivity, colossal mag- 
netoresistance and ferroelectricity [lHj]. Stoichiometric 
SrTiC>3 has an indirect band gap of 3.25 eV and is there- 
fore a band insulator However, its electronic prop- 




; cation-substitutions 
or oxygen vacancies 



erty can be changed by dopin 
like La for Sr or Nb for Ti 

[1, introduce n-type carriers, making STO metal- 
lic. Doping levels that create a carrier concentration as 
low as 2.7 x 10 17 cm~ 3 are sufficient to transform insu- 
lating STO into a metallic state [l5j |. For certain dop- 
ing levels, even superconductivity with T c < 0.3 K has 
been reported [16]. In particular, oxygen vacancies con- 
tribute more than just electric carriers. They exhibit, if 
in great amounts, dark-blue fluorescence through a tran- 
sition between their impurity level and the valence band 
Oxygen vacancies are also actively investigated for 
resistive memory applications [H, [H, Qjj . 

In recent years, researchers found that when thin films 
are grown on STO substrate by pulsed laser deposition 
(PLD), energetic species during the growth process in- 
duce oxygen vacancies and create a metallic layer into 
the substrate [2(|. Similarly, ion-milling process, a com- 
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mon dry-etching technique, was also found to introduce 
metallic states on insulating STO by creating oxygen va- 
cancies [l?], HHHII- These observations raise serious ques- 
tions on transport measurements for films processed on 
STO substrates. Considering the routine usage of ener- 
getic processes such as PLD and ion milling with STO 
substrates, proper understanding of how energetic par- 
ticles create conducting states on STO is important. In 
order to investigate how metallic states emerge on STO 
surface through energetic bombardment, it is essential 
to carry out the experiment without exposing samples 
to air because oxygen in air can quickly fill the vacan- 
cies on the surface and mask metallic states underneath. 
Such studies do not exist yet. Here, by combining ion- 
milling process, in situ conductance measurement and 
model calculations, we show that metallic states emerge 
even when thermal diffusion is completely suppressed at 
cryogenic temperatures, and also show how such conduct- 
ing states emerge and evolve during and after ion milling 
in different conditions. 



II. EXPERIMENTAL SETUP 

The experiment was carried out in a custom-designed 
high vacuum chamber (base pressure less than 10 -7 Torr) 
that is equipped with an in situ four point probe having 
springy contacts; see Ref.[23j for the details of the setup. 
This unique system allows Ar-ion bombardment over a 
wide range of temperatures (— 160 °C to 700 °C) and 
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four-point-probe conductance measurements within the 
same vacuum condition. We used a gridlness ion source 
(EH200F, Kaufman & Robinson), which allows large ion 
beam currents even at very low beam energies down to 
50 eV. We also measured oxygen partial pressures using 
a residual gas analyzer (RGA), and at the base pres- 
sure, the oxygen partial pressure was found to be less 
than 10 -9 Torr. For each experiment, a polished single- 
crystal (001) STO substrate of 10 mm x 10 mm x 0.5 mm 
(MTI and Crystec as received) was mounted with silver 
paint on a copper stage, which was then cooled with ei- 
ther liquid nitrogen or water depending on the target 
temperatures. During ion-milling process, the four point 
probe was positioned far away from the sample so that it 
does not interfere with the ion beam. After each ion- 
milling process, contact between the four point probe 
and the sample was made using a manual linear-motion 
fccdthrough, and the conductance was measured using 
a source meter (Kcithlcy 2636A). Over the entire pro- 
cess, temperature of the sample was monitored using a 
thermocouple mounted on a corner of the substrate. 



III. RESULTS 

A. Model at cryogenic (100 ~ 200K) temperatures 

Because each non-clustered oxygen vacancy donates 
two electrons and the electronic mobility is fairly con- 
stant over a wide ran ge o f doping level at temperatures 
above 100 K in STO [2J], to the first approximation we 
can assume that the electric conductivity of STO mea- 
sured by the four point probe is proportional to the den- 
sity of oxygen vacancies. With this assumption, we first 
construct a simple model on how oxygen vacancies form 
during ion-milling process when thermal diffusion is ne- 
glected and will use it to describe the conductance mea- 
surements at cryogenic temperatures in the following sec- 
tion. 

Oxygen is not the only type of ion that is removed 
during ion milling of the STO surface. In fact, the whole 
surface gets removed, which gives rise to a question why 
oxygen vacancies are formed at all [13, [H| ■ Formation 
of oxygen vacancies can be understood by considering 
the diffusivity difference between ions; the diffusivity of 
oxygen is significantly higher than those of titanium and 
strontium for multiple reasons. First of all, oxygen is 
lighter and is thus easier to move in the lattice. Second, 
its nearest neighbor distance is smaller. The distances 
of Ti and Sr are by a factor of V2 larger than that of 
oxygen. Furthermore, oxygen has a free line of sight to 
the nearest neighbor, which is not the case at least for 
titanium. Because diffusivity depends almost exponen- 
tially on these parameters, their combined effect leads 
to a significant difference in the diffusivities of oxygen 
vs. strontium and titanium. Higher diffusivity of the 
oxygen ion results in the formation of a oxygen vacant 
layer deeper than argon can penetrate into STO. This 
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FIG. 1: (Color online) Development of an oxygen deficient 
layer beyond the penetration depth of argon ions. The sur- 
face layer is continuously removed during argon-ion milling 
process. Oxygen vacancies tend to diffuse deeply into STO 
surface because of the significantly higher diffusivity of oxygen 
than those of strontium and titanium. 



scenario is depicted in Fig. [T] where argon-ion penetra- 
tion depth is smaller than 3 nm for beam energies below 
300 eV [HI, [22| . As the top surface layer is being etched 
away by incoming argon ions, the oxygen ions below that 
layer are able to fill out the vacancies within the penetra- 
tion depth by thermal migration. Through this process, 
oxygen vacancies can be formed in a region much deeper 
than the penetration depth of argon ions [22| . 

Considering that the diffusivity is exponentially depen- 
dent on temperature, the thermal state of the SrTiO,3 
sample during ion-milling must be a critical factor for 
the depth profile of the oxygen vacancies. Based on the 
diffusion coefficients available in the literature for oxygen 
vacancies [25J, at temperatures below — 100 °C, the dif- 
fusion rate for oxygen vacancies is estimated to be much 
smaller than the ion milling rate. This implies that ther- 
mal diffusion can be neglected at these temperatures be- 
cause incoming argon ions will etch away oxygen vacan- 
cies before they diffuse into the bulk. If such is the case, 
from Fig. [T]the oxygen deficient layer should be limited 
to the argon penetration depth, and the density of oxy- 
gen vacancies must reach a steady-state value almost im- 
mediately after ion-milling process is started. However, 
it was not consistent with our experimental observation 
described later, which strongly suggests that the oxygen 
deficient layer extends much beyond the argon penetra- 
tion depth even at these cryogenic temperatures. We 
propose that the lattice vibrations on the surface, caused 
by collisions with the argon ions, are the reason for the 
oxygen vacant layer developing beyond the penetration 
depth of argon ions. The intensity of the vibrations de- 
cays exponentially with depth, and the thickness of the 
vibrating layer, which we call hot-zone, is determined by 
the kinetic energy of the incoming argon ions. Because 
the incoming argon ions are incoherent, coherent vibra- 
tion modes, i.e. phonons, cannot propagate through the 
bulk. The hot-zone formation is also different from the 
standard heating process in that it cannot propagate into 
the crystal beyond the vibration decay length. 

Considering that the propagation rate of oxygen va- 
cancies through the hot-zone must be much higher than 
the ion-milling rate, during short time At, we can as- 
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sume that the oxygen vacancies propagate through the 
hot-zone almost immediately before any significant ion 
milling occurs. If we assume an exponentially decaying 
depth profile of the oxygen vacancies during this short 
time, the depth profile of the carrier concentration, which 
is assumed to be proportional to that of the oxygen va- 
cancies, should increase by 



/ a , dn s A 
An 0, At) = At e 



(1) 



where ^f- is the rate of concentration increase 

at 



at the 

very surface and do is the effective hot-zone thickness. 
Additional factor to consider is the removal of the top- 
most layer during this time. With the milling rate r, the 
thickness of the removed layer should be Ax = r ■ At 
after a time interval At. After each At, new carriers are 
created on the surface (part 2) with the same profile as in 
Eq. (fTJ) and added to the already existing carriers (part 
1) as shown in Fig. [5] If we integrate over the whole 
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FIG. 2: (Color online) Formation of the oxygen vacancy layer 
during ion-milling process, (a) At time t — 0, the oxygen va- 
cancy distribution forms an exponential shape following the 
decay of the vibration energy in the hot-zone, (b) After a 
time interval At, the ion-milling process removes a layer with 
thickness r ■ At, where r is the milling rate. So the vacancy 
concentration profile moves about this depth further in addi- 
tion to the remaining part of the previous profile, (c) As this 
process continues, the original profile (t = 0) gets gradually 
removed. 

depth taking At infinitesimally small, this model leads 
to a areal carrier concentration of 
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Here, the exponential term vanishes after long time and 
the concentration saturates to a steady-state value; note 
that the saturation will occur faster as the milling rate 
r increases. If there were not any etching process (that 
is, r = 0), the areal carrier concentration introduced per 
time, djV "°~7/ c '"" s w0lu d reduce to ^§f-do, which implies 
that the higher the beam current is (larger % 2 -) or the 
higher the ion energy is (larger do), the more vacancies 
are created. It is reasonable to assume that this value 
should be proportional to the milling rate r, because 
milling is nothing more than vacancy formation of all 
the component ions. In other words, 
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no— etching 
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Putting this expression into Eq. ([2]) leads to the simpli- 
fied equation 



N(t) oc d a (l -e~**) 



(4) 



(2) 



This result provides an important insight: the saturation 
carrier density (as t — > oo) is proportional to do, which 
should grow with the argon energy just as the ar gon pen- 
etration depth increases with the beam energy [26[. So 
we expect that higher beam energy will lead to higher 
saturation value of the areal carrier density on ion-milled 
STO surface at cryogenic temperatures. 



B. Measurement at cryogenic temperatures 

Figure [3J shows how electric conductance of fresh in- 
sulating SrTiOa substrates changes as a function of ion- 
milling time with the sample at cryogenic temperatures 
(below — 100 °C), for beam energy of 66 eV and two dif- 
ferent beam current densities. The saturation of the 
conductance after about an hour is consistent with the 
cryogenic-temperature model, Eq. ([2]). According to the 
model, the thickness of the conducting hot-zone, do, must 
be the same for the two curves in Fig. [3J because it de- 
pends only on the beam energy. The fitting parameter 
for the exponent ^ of the I.5mA/cm 2 curve is 2.5 times 
higher than that for the 0.5 m A/cm 2 curve. Despite the 
simplicity of our model, this value is in reasonable agree- 
ment with the factor of three ratio of the beam currents, 
which is also proportional to the ion-milling rate, r. 

However, the fact that the saturation values of the two 
curves are significantly different is inconsistent with Eq. 
(U), according to which the same saturation value is ex- 
pected for both curves because they have the same beam 
energy. We propose that this discrepancy is partly due to 
clustering of oxygen vacancies; this clustering idea is fur- 
ther supported by follow-up measurements discussed be- 
low. When oxygen vacancies form clusters, they start to 
trap electrons, lowering the conductivity from the value 
predicted from non-clustered vacancies [27j . Because 
clustering is energetically favored (28|, unless vacancy 
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FIG. 3: (Color online) Conductance of SrTiOs versus ion- 
milling time with a beam energy of 66 eV and different ion 
beam currents. The ion-milling temperature never exceeded 
— 100 °C, and the conductance measurement was performed 
at — 140 °C. The fitting curves are based on Eq. ([2]). 



thicker oxygen deficient layer, thus higher conductance. 
This is seen in Fig. [U which shows that the conduc- 
tance of the sample can be changed back and forth by 
switching the beam energies. The conductance depends 
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diffusion is completely forbidden, it should be impossi- 
ble to avoid the formation of clusters. Even at cryogenic 
temperatures as in Fig. [3J although thermal diffusion be- 
yond the hot-zone is suppressed, diffusion will still occur 
within the hot-zone during ion milling, and thus nearby 
vacancies can form clusters. The longer vacancies stay 
within the hot-zone during ion milling, the more clusters 
will form, resulting in lower electrical conductivity. This 
scenario is consistent with Fig. [3J The crucial factor is 
the time it takes to refresh the hot-zone. This time can be 
calculated out of the milling rate r and the thickness do- 
The faster milling rate of the 1.5mA/cm 2 measurement 
leads to a higher refreshing rate, so the oxygen vacancies 
have less time to cluster and the electric conductance 
should be higher than that of the 0.5mA/cm 2 measure- 
ment. 

The thickness of the conducting layer, say for the 
I.5mA/cm 2 curve, is estimated as do — 57 nm from the 
fitting parameter and measured milling rate. It should 
be noted that this thickness is that of the oxygen defi- 
cient layer and is much larger than the ar gon penetration 
depth, which is only a few nanometers [21| . Based on 
the effective thickness of the conducting layer, do, and 
the electron mobility (/i = 49 ^-) of conducting STO at 
the measurement temperature (—140 °C) [29|, the steady- 
state surface carrier density, n s (= is estimated 
as n s .66eV — 2.0 x 10 20 cm -3 , which is comparable to high 
end values due to oxygen vacancies [29j | . This is equiva- 
lent to about 0.3% oxygen vacancies at the top surface, 
if every oxygen donates two electrons. 

According to our model, if the beam energy is changed, 
thickness of the hot-zone, accordingly the conductance, 
should also change; higher beam energies should result in 



FIG. 4: (Color online) Conductance switching of SrTiOs at 
— 130 °C with alternating ion beam energy and constant beam 
current density of 0.5m A/cm 2 . Higher energy provides a 
thicker hot-zone, and thus higher conductance. Switching 
back to 66 eV causes the conductance to decrease as the con- 
ducting layer is getting thinner. Such switching behavior did 
not occur near room temperature. 

on the beam energy because the thickness of the hot-zone 
changes. It is noted that the observed behavior is a com- 
bination of various effects. The influence of clusters and 
increasing surface roughness are factors that complicate 
the development of an exact model, but it is intriguing 
that our simple model can still explain the main features 
of the observed phenomena. 



C. Model near room temperatue 

Measurements near room temperature showed a totally 
different behavior. It was found that the thermal dif- 
fusion of oxygen vacancies can be described by Fick's 
second law [251 ] . The significantly enhanced mobility of 
vacancies at room temperature compared to that at cryo- 
genic temperatures leads to much thicker conducting lay- 
ers with reduced carrier densities. Nevertheless, it is still 
reasonable to assume that a steady-state carrier concen- 
tration, n s , should be reached at the top surface at a 
certain time after ion milling initiates; this will be the 
case when the amount of vacancies migrating from the 
surface into the bulk equilibrates with the amount of va- 
cancies introduced onto the surface by ion-milling. This 
constraint of a constant surface concentration leads to 
the following time and space dependent solution of Fick's 
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second law 



n (x, t) = n s erfc 



(— 

\2Dy/t 



(5) 



where D is the vacancy diffusion constant, x is the depth 
from the surface, and erfc is the complementary error 
function. This equation describes development of the 
carrier density per volume with increasing depth x and a 
constant surface concentration n s . Integrating this func- 
tion over the entire thickness leads to the following areal 
carrier concentration: 



JVtotal (*) 



2n s y/~D 



Vt 



(6) 



the 200 eV. The estimated surface concentrations are 
n s ,66cv — 1-8 x 10 18 cm~ 3 and n Si2 oocV = 1-5 x 10 18 cm -3 
for the 66 eV and the 200 eV case, respectively, which 
are more than two orders of magnitude smaller than the 
value, n S fi &e v = 2.0 x 10 20 cm~ 3 , found above at cryo- 
genic temperatures. These numbers are all physically 
quite reasonable and suggest that both our low temper- 
ature model with hot zone and high temperature model 
with thermal diffusion explain the main mechanisms be- 
hind ion-milling induced metal-insulator transitions on 
STO. 



E. Thermal cycle and Time dependence 



It is noted that this areal carrier concentration never sat- 
urates, unlike the cryogenic case above; this implies that 
the oxygen vacancies contiously diffuse into the bulk at 
these elevated temperatures. 



D. Measurement near room temperatue 

The experimental results are consistent with this dif- 
fusion model as shown in Fig. [5] The plot shows the 
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FIG. 5: (Color online) Conductance versus milling time for 
two different beam energies near room temperature. The be- 
havior is completely dominated by thermal diffusion. The fit 
is done with Eq. [6] The slightly higher conductance of the 
200 eV curve is accounted for by the higher sample temper- 
ature during ion milling induced by its higher beam energy. 
The conductance measurement was done at 30 °C. 

development of the conductance with different beam en- 
ergies and with the same beam current. The fit with Eq. 
|5] is valid only after a certain time because it takes time 
to reach the steady-state surface concentration, which is 
the assumption used for Eq. [6] The effective thickness 
of the oxygen deficient layer increases over time with 
Vt, and, at t = lOOmin, is estimated as d c ff,66cV = 
2.4 /zm for the 66 eV case, and d c ff,200ev = 3.7 /zm for 



Figure [5] shows how the conductance of a cryogeni- 
cally ion-milled sample changes over a thermal cycle up to 
room temperature. It shows that an irreversible change 
of carrier concentration occurred during the cycle, and 
this phenomena cannot be explained by our simple mod- 
els assuming two electrons per oxygen vacancy. We pro- 
pose that vacancy clustering should be the most likely 
mechanism behind this effect. In Fig. [SJ a fresh in- 
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FIG. 6: (Color online) Conductance change due to thermal 
cycle up to room temperature. After ion-milling at cryo- 
genic temperatures, the sample conductance was contiually 
measured during the warm-up and cool-down. Note that 
the low temperature conductance significantly (factor of four) 
dropped after the thermal cycle. In contrast, maintaining the 
sample at cryogenic temperatures preserved the conductance 
for over 20 hours. 

sulating STO was ion-milled at a cryogenic tempera- 
ture (below — 100 °C) and its conductance was measured 
through a thermal cycle from — 180°C up to room tem- 
perature. As the sample is warmed to room temperature, 
its conductance should drop because the electron mobil- 
ity decreases in STO as temperature rises [2{|. If this is 
the only mechanism affecting the conductance, the ini- 
tial conductance should be recovered when the sample 
is cooled back to its initial temperature. Instead, the 
conductance dropped by a factor of four. 
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There are two mechanisms that can explain this irre- 
versible process: one is oxygen vacancy clustering and the 
other is surface oxidation. Both are diffusion-limited, and 
thus suppressed at cryogenic temperatures, but as the 
sample warms up, both mechanisms become active and 
cause the conductance to decrease. Considering that the 
entire experiment was performed in a high vacuum cham- 
ber with oxygen partial pressure less than 10~ 9 Torr, the 
oxidation effect is less likely and thus clustering seems to 
be the dominant effect causing this phenomenon. Ther- 
mal diffusion, which is negligible at cryogenic temper- 
atures, enables oxygen vacancies to migrate and form 
clusters near room temperature. This state cannot be re- 
versed upon re-cooling because clusters are energetically 
favorable [HI, [3(1 EiH . When two vacancies join to form 
a cluster, only half of the original carriers remain mo- 
bile due to the trapping mechanism 28]. Accordingly, 
clustering with two vacancies per cluster would yield a 
conductance drop of only 50% even if the entire vacan- 
cies have clustered. The observed factor-of-four conduc- 
tance drop, therefore, implies that significant portion of 
the clusters were composed of more than two vacancies 
per cluster. 

Even though the partial oxygen pressure during the 
measurement was less than 10~ 9 Torr, it may still be pos- 
sible for oxygen vacancies to be filled by oxygen species 
from the vacuum chamber. If such oxidation occurs, the 
conductance would drop even without the clustering ef- 
fect. In order to distinguish the two effects, we carried 
out an experiment shown in Fig. [7] It shows how con- 
ductance changes in an oxygen-deficient STO surface as 
the sample is exposed to increasing oxygen pressure over 
a long period of time. Initially, the conductance showed 
a fast drop and then it changed to a slow decrease after 
about five hours. During that time, the partial oxygen 
pressure remained below 10~ 9 Torr. By increasing the 
oxygen pressure gradually to 510 Torr at around 50 hours, 
the conductance showed another fast drop followed by the 
same slow decrease as before. This slow background be- 
havior did not change even after pumping out the oxygen 
down to 10 -9 Torr level after 120 hours. This clearly im- 
plies that there are at least two different mechanisms that 
can reduce the conductance. The sudden conductance 
drop, seen after introducing oxygen at around 50 hours, 
is clearly due to oxygen filling up the vacancies. How- 
ever, the overall slow conductance decrease regardless of 
the background oxygen pressure cannot be explained by 
oxidation, and the only known mechanism that can ex- 
plain such a phenomenon is the vacancy clustering. 

Because clusters are energetically favorable @, [28| , dif- 
fusion will force oxygen vacancies to migrate to form 
clusters, localizing otherwise itinerant electrons around 
them and causing the conductance to drop. The amount 
of mobile single vacancies will decrease as the clustering 
develops, and the clustering rate, thus conductance drop, 
will be high at the beginning and will reduce as time goes 
by. Therefore, the initial fast drop before 5 hours in Fig. 
[7] could also be explained by clustering even if oxygen fill- 
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FIG. 7: (Color online) Conductance decrease of an ion-milled 
STO surface over a long period of time with exposure to in- 
creasing oxygen pressure. Both the ion-milling and the mea- 
surement were carried out at room temperature. Note that 
the gradual conductance decrease before and after 50 hours 
and after 120 hours remain almost unaffected even with orders 
of magnitude change in the oxygen pressure. 



ing of the vacancies were completely forbidden. Although 
it is hard to separate out the vacancy-filling contribution 
to the time-dependent conductance drop, the important 
message is that vacancy clustering, which reduces the 
amount of mobile carriers for each oxygen vacancy, oc- 
curs actively in STO even at room temperatures. 



IV. CONCLUSION 

In conclusion, the ion-milling process, a common dry- 
etching technique for complex oxide samples, creates on 
SrTiOs surface a highly conducting layer much deeper 
than argon penerates, not only near room temperature 
but also at cryogenic temperatures where thermal dif- 
fusion is completely suppressed. At cryogenic tempera- 
tures, the conductance gradually grew and saturated to 
a highly metallic value, as a steady-state oxygen defi- 
cient layer developed within hot zone created by bom- 
bardment of the argon ions. Near room temperature, 
however, the condutance kept increasing without satu- 
ration because themal diffusion allowed oxygen vacan- 
cies to diffuse continuously into the bulk. These metal- 
insulator transitions driven by ion-bombardment process 
were well described by simple models assuming each oxy- 
gen vacancy donating two mobile electron carriers. The 
discrepancy between the measurements and the models 
can be accounted for by vacancy clustering, which makes 
oxygen vacancies less effective in donating electric carri- 
ers. Whether ion-bombardment-driven metal-insulator- 
transitions exist in materials other than SrTiC>3 or not is 
an interesting open question for future studies. 
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